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I.   INTRODUCTION 

The  Sylvania  Model  941  CO2  Laser  to  be  used  in  atmospheric 
extinction  studies  was  determined  to  have  poor  beam  quality 
and  frequency  stability  of  the  laser  output.   Since  the 
chance  of  oscillation  on  a  specific  vibrational-rotational 
line  is  dependent  on  whether  a  cavity  mode  is  in  resonance 
with  a  vibrational-rotational  peak,  the  mechanism  used 
in  frequency  tuning  a  CO_  laser  involves  varying  the  cavity 
mode  spacing  by  changing  the  cavity  length.   As  the  laser  is 
tuned  in  frequency  by  changing  the  mirror  separation, 
competition  effects  occur  and  oscillation  can  switch  from 
one  rotational  level  to  another.   The  fitting  of  a  piezo- 
electric driver  to  one  of  the  laser  cavity  mirrors  has  made 
possible  the  use  of  active  feedback  control  to  frequency 
stabilize  the  output. 

The  use  of  active  stabilization  to  minimize  frequency 
drift  has  been  approached  in  a  number  of  ways.   It  requires 
a  comparison  of  the  laser  output  frequency  with  a  suitable 
frequency  standard,  generating  a  correction  signal  to  be 
fed  back  to  the  piezoelectric  transducer  to  adjust  the  laser 
cavity  optical  length.   Such  a  signal  may  be  obtained  by 
passing  a  fraction  of  the  laser  output  intensity  to  a 
suitable  detector  through  an  absorption  cell  containing 
a  suitable  absorbing  gas.   The  Stark  effect  permits  tuning 
an  absorption  line  of  NH2D  to  oscillate  about  the  desired 
laser  vibronic  transition  frequency.   The  detector  then 
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gives  a  sign-dependent  error  signal.   The  ability  to 
generate  a  discriminant  voltage  without  impressing  any 
modulation  on  the  output  and  the  voltage  control  of  the 
frequency  reference  make  the  Stark  cell  stabilization 
technique  attractive  as  a  control  system. 

This  thesis  describes  the  continuation  of  studies  of 
the  operation  of  Stark  cell  modulation,  and  the  technical 
requirements  for  stabilization  of  a  carbon  dioxide  laser. 
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II.   THEORETICAL 

A.   CARBON  DIOXIDE  LASER 

1.   Lasinq  States  and  Spectra  of  Carbon  Dioxide 

The  quantum  states  and  energy  levels  on  which 
molecular  lasers  operate  originate  from  three  contributions: 
i)  electronic  energy,  due  to  the  electronic  charge  cloud 
surrounding  the  nuclei;   ii)  vibrational  energy,  due  to  the 
motion  (vibrations)  of  the  nuclei  about  their  equilibrium 
positions;   and  iii)  rotational  energy,  due  to  the  rotational 
motion  of  the  molecule  about  various  axes  in  space.   Each 
of  these  components  is  quantized.   Each  widely  spaced 
electronic  level  is  split  into  a  series  of  more  closely 
spaced  vibrational  levels.   This  energy  level  scheme  is 
shown  in  Figure  1.   Molecular  gas  lasers  exploit  transitions 
between  the  vibrational-rotational  levels  of  a  molecule. 

In  general,  an  N-particle  system  has  3N  degrees  of 
freedom.   The  motion  of  the  center  of  gravity  and  the 
orientation  of  the  axis  of  a  linear  XY  molecule  can  be 
specified  in  terms  of  five  degrees  of  freedom,  leaving 
(3N  -  5  =  4)  four  vibrational  degrees  of  freedom.   For  a 
linear  triatomic  symmetric  molecule  such  as  CO  ,  there  are 
three  distinct  normal  modes  of  vibration.   The  four  vibrations 
correspond  to  the  symmetric  stretch  mode,  the  doubly  degen- 
erate bending  mode  and  the  as^^metric  stretching  mode.   These 
modes  are  shown  in  Figure  2 .   In  a)  the  symmetric  stretch 
mode,  the  oxygen  atoms  vibrate  along  the  internuclear  axis 
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in  a  symmetric  manner.   In  b)  the  bending  mode,  the  atoms 
vibrate  symmetrically  in  a  direction  perpendicular  to  the 
internuclear  axis.   The  bending  mode  has  a  two- fold  degeneracy 
which  arises  from  the  equivalence  of  vibration  in  either  of 
two  orthogonal  planes  whose  intersection  lies  along  the 
molecular  axis.   In  c)  the  asymmetric  mode,  the  atoms  vibrate 
asymmetrically  along  the  internuclear  axis.   Since  carbon 
dioxide  is  symmetrical  in  its  equilibrium  state  it  does  not 
have  a  dipole  moment.   The  symmetrical  stretching  vibration 
does  not  create  a  net  electric  dipole  moment  and  therefore 
is  not  active  in  the  infrared.   The  asymmetric  normal 
vibrations  produce  a  dipole  moment  and  are  I.R.  active. 
The  stretching  vibrations  have  higher  frequencies  than  the 
bending  vibrations  since  the  restoring  force  is  greater  in 
stretching  a  molecule  than  in  bending. 

The  C0_  molecule  vibrational  quantum  state  can  be 
described  by  a  set  of  quantum  numbers  written  as 


[\-  "r  "3] 


where  v  is  the  degree  of  excitation  of  each  of  the  three 
individual  normal  modes .   The  degenerate  bending  vibration 
mode  for  CO  has  an  associated  angular  momentum  which 
depends  upon  the  excitation  and  relative  phase  of  the 
orthogonal  vibrations.   This  angular  momentum  is  denoted 
by  the  quantum  number  JZ  appearing  above  v  .   The  superscript 
i     takes  the  value  i    =   v  ,  v   -2,...,  1,  0,  with  each 
giving  rise  to  a  sublevel. 


Assuming  independence  of  the  three  normal  modes,  and 
the  use  of  normal  coordinates  in  describing  the  molecular 
vibrations,  allows  the  separation  of  the  Schrodinger  equation 
into  as  many  separate  equations  as  there  are  normal  modes. 
The  eigenvalues  for  each  of  the  separate  equations  are 
given  by 

Ei  =  (v^  +  h)    hcVj_  (1) 


where  v.  is  the  quantum  number  associated  with  the  i*-"  mode, 
and  the  total  vibrational  energy  associated  with  an 


rl'  ^2'  ^3  1  st^"^®  ^s  approximately 


E^  =  (v-L  +  h)    hcV^  +  (V2  +  h)    hcV2  +  (V3  +  h)    hcV3     (2) 

where  v  is  the  frequency  in  wave  number  units  associated 
with  the  normal  modes  of  vibration. 

The  carbon  dioxide  vibrational  energy  levels  of 
principal  interest  in  laser  action  with  their  modes  and 
frequencies  are 

1)  [oo^l"]  Asymmetric  Stretch,  2^  2  349.3  cm 

2)  [lO°OJ  Symmetric   Stretch  1333.3    cm"-'" 

3)  [o2°oJ  Bending, Z!g  1235.5    cm' 

4)  roi-'"o']  Bending  667.3   cm 


-1 
-1 
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The  two  vibrational  levels  [lO°oJ  and  |_02°oJ  belonging  to 
different  vibrational  modes  have  nearly  the  same  energy 
and  are  in  Fermi  resonance.   When  the  interaction  between 
them  is  included  as  a  perturbation,  the  wave  functions  are 
mixed  and  the  resulting  pair  of  mixed  levels  are  repelled, 
so  that  the  level  separation  is  greater  than  expected. 

The  lower  vibrational  energy  level  diagram  for  the 
CO2  molecule  in  conjunction  with  the  N2  molecule  is  shown 
in  Figure  3.   The  dashed  lines  indicate  radiative  transitions 
that  are  allowed  by  electric  dipole  considerations.  The  Zl 
and  TT  notation  represents  the  electronic-vibrational 
angular  momentum  of  the  molecule  (2=0,'n"=i),  giving  the 
value  of  the  vector  sum  of  the  electronic  orbital  angular 
momentum  and  vibrational  angular  momentum  about  the  inter- 
nuclear  axis .   The  +  superscript  implies  that  the  vibronic 
eigenfunction  is  unchanged  when  reflected  at  any  plane 
containing  all  the  nuclei.   the  subscripts  indicate  whether 
the  state  is  u,  antisymmetric  (odd  symmetry)  or  g,  symmetric 
(even  symmetry)  with  respect  to  vibrational  motion.  jjRet .  i] 

For  every  vibrational  state  there  exists  a  distinct 
set  of  rotational  levels  with  slightly  different  spacing 
for  different  vibrational  levels.   The  CO2  molecule  can 
rotate  about  any  axis  normal  to  the  internuclear  axis. 
From  quantum  theory  of  rotational  spectra  and  angular  momentum 
the  energy  level  due  to  rotational  motion  is  given  by 

E_(J)  =  BJ(J  +  1)  =  h J(J  +  1)        (3) 
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where  J  is  the  angular  momentum  quantum  number  (J  =  1,  2, ...) , 

I  is  the  moment  of  inertia,  and  B  is  called  the  rotational 

constant.   The  selection  rule  for  J  for  dipole  transitions 

of  COp  allows  only  transitions  between  rotational  levels 

for  which A  J  =  -1. 

Laser  action  for  carbon  dioxide  occurs  on  transitions 

between  two  vibrational  levels  of  the  same  electronic  state 

or  more  precisely,  transitions  between  rotational  levels  of 

the  two  vibrational  states  (vibration-rotational  transitions) . 

These  transitions  are  classified  into  "branches"  by  the 

change  in  J.   Transitions  between  two  vibrational  levels  v 

and  V  with  account  taken  for  the  rotational  splitting  are 

shown  in  Figure  4.   The  transition,  in  the  absence  of 

rotational  energy,  would  consist  of  a  single  line  centered 

at  \)    ,    called  the  Q  branch  where  AJ  =  0.   Including  rotational 
^o 

energy  the  transition  actually  consists  of  two  groups  of  lines: 
the  P  branch  corresponding  to  a  change  in  rotational  quantum 
number  of  AJ  =  +1  and  R  branch  corresponding  to  a  change  in 
rotational  quantum  number  of  A J  =  -1. 
Solving  for  V  in  the  equation 

AE  =  hV  (4) 


and  defining  v   ^s  the  band  center  frequency  or 
o 


Vo  =  F  K    '°>  -  V  '°il  <5) 
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the  transition  frequencies  can  be  written  as 


V  =  V   -(^^/h)J  P(J  +  1)  Branch 

o 

Vq  =  V^  Q  Branch  (b) 

V  =  V^  +  (^^/h) J  R(J  -  1)  Branch 


where  J  is  the  rotational  quantum  number  of  the  lower  level. 
The  rotational  levels  of  a  linear  molecule  are  positive  or 
negative  depending  on  whether  the  sign  of  the  total  eigenfunction 


^r='i'e^r'^v  '"' 


remains  unchanged  or  changed  upon  an  inversion. 

The  infrared  selection  rules  for  the  vibrational- 
rotational  transitions  relevant  to  CO  are  [kef.  2~1 


Av  =  -1 


A.i  =  0,  -1;   g-z-g;   u -^Z*  u  (3) 

A  J  =  0,  tl;   +-.-*.+;   s-^V-  a   (J  =  0)-»7^  (J  =  0) 


where  s  and  a  designate  the  symmetric  and  antisymmetric 

total  rotational  eigenfunction  "iL     and  the  symbols  -•— ►  and 

■^r 


represent  allowed  and  not  allowed  respectively.   In  all 
symmetric  excited  vibrational  levels  (Zl  )  the  even  rotational 
levels  are  positive  and  the  odd  are  negative  with  respect  to 
inversion  along  the  molecular  axis  while  the  converse  is  true 
for  the  antisymmetric  levels  (Zi~).   For  the  TT  ,  A  ...  doubly 
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degenerate  vibrational  levels,  for  each  J  there  is  a  positive 
and  negative  level  with  alternating  order  (+-,  -+,  +-,  -+,... 
or  -+,  +-,  -+,  +-,...). 

Because  CO^  is  exactly  symmetric  and  is  without 
nuclear  spin  the  Q  branch  or  AJ  =  0  transition  is  strictly 
forbidden  by  symmetry  restrictions  (ZT*—*  Z!  /  ^J  7^  0)  .   J  to 
J  ;tl  transitions  within  a  single  vibrational  energy  level 
are  also  forbidden.   Since  the  two  oxygen  atoms  are  spinless, 
symmetrically  located  indistinguishable  particles  and  the 
CO2  wave  function  remains  unchanged,  the  system  obeys  Bose- 
Einstein  statistics.   Therefore  all  CO  g  states  have 
rotational  levels  with  even  J  values  only,  while  all  u  states 
have  only  odd  J  values.   The  absence  of  alternating  rotational 
levels  causes  the  spacing  between  adjacent  lines  to  be  twice 
that  expected  on  the  basis  of  the  size  of  the  molecular 
rotational  constant  B.   This  spacing  is  4B  instead  of  2B. 
Experimental  evidence  shows  that  only  the  even  J  values 
appear  in  CO   P-  and  R-branch  spectra.   The   C  O2  system 
can  oscillate  on  approximately  one  hundred  transitions 

spanning  the  wavelength  range  from  9 . 1  to  lO.l^m-      They 

o       o 
occur  on  the  P(J)  and  R(J)  branches  of  both  the  00  1*—*  10  0 

and  00  1"*-*  02  0  rotational-vibrational  bands.   A  detailed 

o      o        00 
laser  transition  diagram  for  the  00  1  -  10  0  and  00  1  -  02  0 

bands  including  some  rotational  levels  is  shown  in  Figure  5 

o 
and  some  of  the  P-  and  R-branch  transitions  for  the  00  1  - 

10  0  vibration-rotation  levels  with  their  corresponding 

transition  wavelengths  and  frequencies  are  given  in  Table  I. 

[Ref.  3] 
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TABLE  I 


00  1  -  10°0  TRANSITION 


TRANSITION 
P(14) 
P(16) 
P(18) 
P(20) 
P(22) 
P(24) 
P(26) 
P(28) 


P  BRANCH 
WAVELENGTH  (  m) 
10.5326 
10.5513 
10.5713 
10.5912 
10.6118 
10.6324 
10.6534 
10.6748 


FREQUENCY  (cm~  ) 
949.43 
947.70 
945.96 
944.18 
942.35 
940.52 
938.67 
936.78 


TRANSITION 
R(14) 
R(16) 
R(18) 
R(20) 
R(22) 
R(24) 
R(26) 
R(28) 


R  BRANCH 
WAVELENGTH  (  m) 
10.2360 
10.2855 
10.2605 
10.2470 
10.2335 
10.2200 
10.2075 
10.1950 


FREQUENCY  (cm~  ) 
971.91 
972.24 
974.61 
975.90 
977.18 
978.47 
979.67 
980.87 
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2 .   Excitation  Mechanism 

Electrical  pumping  in  a  gas  laser  is  achieved  by 
allowing  a  current  to  pass  through  the  gas  producing  ions 
and  free  electrons .   The  average  electron  energy  and  the 
ability  of  an  electron  to  excite  a  neutral  molecule  by 
collision  is  much  greater  than  the  ion  energy.   The  equili- 
brium condition  established  among  the  electrons  can  be 

described  by  an  effective  electron  temperature  T  .   Assuming 

e 

a  Maxwellian  energy  distribution,  and  that  for  each  collision 
some  fraction  6  of  the  electron  energy  is  lost,  the  effective 
electron  temperature  is  given  by  [kef.  4j 


T«  =  -7-t4 (Eji  )  (9) 

e     {2  6  f^   k 


where  E  is  the  applied  electric  field,  k  is  Boltzmann's 
constant,  e  is  the  elementary  charge  on  an  electron,  and  i 
is  the  electron  mean  free  path.   Since  the  mean  free  path  is 

i   =  ,-^^^  1 —  (10) 

2 
with  TTo*   the  collision  cross  section  and  P  the  pressure, 

E  / 
T  can  be  given  as  a  function  of  the  ratio   /P.   The  ratio 
e        ^ 

of  applied  electric  field  to  pressure  is  used  as  a  parameter 
in  describing  the  discharge  condition. 

The  primary  electrical  pumping  mechanisms  for  a 
gas  are:   i)  excitation  by  electron  impact  or  a  collision 
of  the  first  king  represented  by  the  equation 
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A  +  (e  +  KE)  — ►   A*  +  e  (11) 

where  *  represents  the  excited  state  of  A,  and  KE  the 
translational  kinetic  energy  converted  to  excitation 
energy.   ii)  Excitation  by  collision  between  atoms  or 
molecules  of  different  species  (resonant  energy  transfer) 
characterized  by  the  equation 

A  +  B*  ►   A*  +  B  +  KE  (12) 

where  A*  and  B*  are  the  excited  states.   This  type  of 
excitation  process  is  called  a  collision  of  the  second  kind. 
Selective  excitation  of  polyatomic  molecules  through 
vibrational  energy  transfer  from  a  diatomic  molecule  has 
been  used  effectively  in  the  production  of  population 
inversion  on  vibrational-rotational  transitions  of  the 
polyatomic  molecule.   For  the  N  -  CO^  system  shown  in 
Figure  3,  the  [OO  ij  level  of  carbon  dioxide  is  chiefly 
excited  by  a  resonant  energy  transfer  via  collisions  of  the 
second  kind  between  CO2  molecules  in  the  ground  state  and 
N_  molecules  in  the  excited  state.   As  a  result  of  these 
resonant  collisions  the  populations  of  certain  sets  of 
levels  are  very  strongly  coupled  so  that  they  act  essentially 
as  a  single  level.   The  vibrationally  excited  N  molecules 
in  the  ground  electronic  state  N  *  (v  =  1)  are  produced 
primarily  by  indirect  electron  impact  excitation.   The 
pump  rate  R  of  the  upper  state  is 


21 


Rd  =   N^      <V(r>  (13) 


T        "e 


where 


<V<r>       =    J*  VO- (E)  f  (E)dE  (14) 


N  is  the  electron  density,  V  is  the  electron  velocity, 

f(E)  is  the  electron  energy  distribution  and  0"  (E)  is  electron 

excitation  cross  section.  [_Ref.  sj   The  cross  section  for 

electron  impact  excitation  of  N2  in  the  ground  state  is 

very  high  making  substantial  population  densities  of  N  *  (v  =1) 

achievable.   Figure  6  shows  the  total  electron  excitation 

cross  section  of  the  N   levels  as  a  function  of  electron 

2 

energy  with  the  cross  section  reaching  a  maximum  value  of 

—  1  fi   9 
approximately  3  X  10   cm  at  an  electron  energy  of  2.3ev. 

[Ref.  6]   The  vibrationally  excited  nitrogen  molecule  in 

the  ground  electronic  state  cannot  decay  through  electric 

dipole  radiation  back  to  the  ground  state  (v  =  1  to  v  =  0) 

since  N„  is  a  diatomic  homonuclear  molecule  with  no  permanent 

dipole  moment.   The  energy  difference   E  between  the  N  *  (v  =  1) 

and  CO^  [j30°l]  energy  levels  is  about  18cm   .   This  is 

negligible  in  comparison  with  the  average  thermal  energy 

kT  at  room  temperature  of  210cm   .   The  near  resonance 

molecule  -  molecule  collision  leading  to  energy  transfer 

can  be  written  as 


N*    (v=l)    +   CO       [OO^OJ N2     (v=0)    +   CO       [00°l]-18cra  (15) 
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This  resonance  transfer  produces  a  mixed  state  in  which  the 
cornbined  populations  of  N  *  (v  =  1)  and  CO^  [00°l]  are 
essentially  in  equilibrium:   this  effect  increases  the 
effective  lifetime  of  CO   [OO  IJ  .   The  transfer  cross 
section  is  largely  due  to  the  negligible  difference  between 
the  energy  levels.   The  level  flO^o]  (or  [p2  o]  )  is  more 
than  900cm   from  N  *  (v  =  1)  and  the  transfer  cross  section 
is  small  compared  to  that  for  the  reaction  in  (15) . 

To  a  lesser  degree  excitation  can  occur  by  direct 
electron  impact  on  the  ground  state  CO2  molecule  [oo'^OJ 
and  by  resonant  exchange  of  energy  between  CO  molecules 
in  excited  levels  and  C0_  in  the  ground  state.   Electron 
impact  on  the  ground  state  CO2  molecules  results  in  excitation 
of  the  CO2  molecule  to  the  [00°lj  level  since  the  transition 
to  ZJ    [00  iJ  is  optically  connected  to  the  ground  state 
Zk         Loo  OJ  and  is  given  by 

CO2  \j^      .    00°0]  +  (e  +  XE) C02[2  ^ ,    00  l]  +  e    (16) 


The  cross  section  for  vibrational  excitation  of  CO  by 

electron  impact  is  shown  in  Figure  7.  [Ref .  vj 

A  substantial  amount  of  CO  (10%  of  CO2  pressure) 

can  be  produced  as  a  result  of  dissociation  of  CO   in  a 

2 

discharge.   The  electron  excitation  cross  sections  for  CO 
are  reported  in  Reference  8  to  be  similar  in  form  to  that 

of  Np  but  shifted  to  lower  energies  with  the  peak  value  for 

-16   2 
the  total  cross  section  being  approximately  3  X  10    cm 
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at  1.8  ev.   The  difference  between  the  first  CO  vibrational 
energy  level  and  the  upper  laser  level  of  CO   is  given  as 
ITOcm   .   These  final  two  mechanisms  account  for  the 
possibility  of  producing  a  population  inversion  in  a  pure 
CO2  laser. 

Normally  the  CO   oscillator  functions  with  a  plasma 
that  consists  of  a  CO^  -  N  -  He  mixture.   The  helium 
component  plays  two  important  roles.   The  first  is  related 
to  the  increased  rate  of  relaxation  of  the  [oi  oj  level. 
By  reducing  the  bending  mode  population  which  in  turn  reflects 
itself  in  a  reduction  of  the  lower  laser  level  populations 
the  gain  is  increased.   Helium  also  reduces  the  kinetic 
and  rotational  temperature  of  the  CO  in  the  plasma.   This 
will  also  increase  the  gain,  since  the  same  total  upper- 
state  population  is  now  distributed  among  fewer  levels. 
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3 .   Lineshape  Function 

The  function  describing  the  distribution  of  emitted 
intensity  versus  the  frequency, V  ,  is  called  the  lineshape 
function,  g(V)  and  is  normalized  according  to 


j^g(V)  dV  =  1  (17) 


-oo 


The  frequency  spectrum  and  g(V)  are  the  same  for 
purposes  of  spontaneous  emission,  stimulated  emission  and 
absorption.   Line  broadening  mechanisms  may  be  divided  into 
two  classifications:   i)  homogeneous  and  ii)  inhomogeneous . 

Homogeneous  broadening  involves  the  spectral  spread 
of  the  whole  system  caused  by  the  broadening  of  the  spectrum 
of  each  indistinguishable  molecule-   It  is  due  most  often 
to  the  finite  collisional  lifetime  of  excited  molecules 
and  the  consequent  interruption  of  the  radiative  process. 
The  result  of  frequent  collisions  is  the  presence  of  many 
truncated  radiative  or  absorptive  processes  with  a  line- 
width  of  radiation  greater  than  that  of  an  individual 
uninterrupted  process. 

The  normalized  Lorentzian  lineshape  function  which 
is  characteristic  of  homogeneous  broadening  is  [Ref.  9] 


3<^^=  airttv-v^t'^^^v./:,)^]  '^'' 


Here  v^^  is  the  center  frequency  and  AV|_  is  the  separation  between 
the  two  frequencies  at  which  the  Lorentzian  is  down  to  half 

26 


its  peak  value,  or  the  full  width  at  half  maximum. 

The  linewidth  AVj_  is  related  to  the  finite  lifetime  T 
of  the  emitting  state  by 

The  Lorentzian  linewidth  is  pressure  dependent 
since  the  frequency  of  collisions  is  proportional  to  the 
density  of  the  gas.   The  time  interval  between  collisions 

can  be  approximated  by  the  ratio  of  the  mean  free  path 
(10)  to  the  mean  thermal  velocity  V 


-    ,      akT   h 

^   =  (  -WW   )  (20) 


and 


(MkT)  2 

^    -ir — -  (21) 


c     {I6ir)'2  Pa 

where  M  is  the  mass  of  the  molecule.   Since  'T.  is  inversely 
proportional  to  the  pressure  P,  the  linewidth  ^V  is 
proportional  to  the  pressure. 

Inhomogeneous  broadening  is  due  to  the  spread  in 
the  individual  transition  frequencies  of  the  molecules  and 
not  the  broadening  of  the  lines  of  the  individual  molecules. 
The  thermal  motion  of  the  gaseous  molecules  is  the  cause 
of  Doppler  broadening.   The  transition  frequency  V  is 
Doppler  shifted  due  to  the  finite  veldcity  of  the  molecule. 
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Using  the  Maxwell  velocity  distribution  function  of  a  gas, 
the  normalized  Doppler  broadened  lineshape  function  for 
inhomogeneous  broadening  is  Gaussian  and  given  in  Reference  10 
as 


9(^'  =  ^  7WW    ""   e-p[-("/2'^T)  i^%l)i^-^y]  (22) 


where  M  is  the  molecular  mass  at  equilibrium  at  temperature  T, 
and  Vq  is  the  resonant  frequency. 

The  linewidth  obtained  from  (22)  is 


AV,=  2V  L^'^-'^Vmc^]'^  (23) 


and  reexpressing  g(V)  in  terms  of  AV_^  gives 


g  (V )  =  1 Exp  I  ° '  /av_  I 


(24) 


For  a  carbon  dioxide  laser  Doppler  broadening 
dominates  at  low  pressure  giving  a  Gaussian  line,  and 
pressure  (collisional)  broadening  dominates  at  high  pressures 
giving  a  Lorentzian  line.   In  general,  the  overall  broad- 
ening is  due  to  a  combination  of  effects  and  mechanisms  and 
is  given  by  a  convolution  of  the  various  processes.   The 
overall  lineshape  (g  )  of  two  broadening  mechanisms,  one 
homogeneous  (g  )  and  the  other  inhomogeneous  (g  ) ,  may  be 
written  as 

oo 

^o  ^^^  "  J  gi(^-V')g2(V')  dv'  (25) 


-oo 
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The  convolution  of  a  single  Lorentzian  line  with 
a  single  Gaussian  in  (25)  gives  a  Voigt  profile.  [Ref.  Ill 
If  both  Doppler  and  collision  widths  are  appreciable,  the 
total  line  width  can  be  approximated  by  [kef.  lij 


AV  '^  [aV^        +^V^  1  ^        (26) 


Figure  8  (a)  and  8  (b)  shows  a  general  comparison 
of  Lorentzian,  Gaussian,  and  Voigt  lineshapes .   In  8  (b) 
the  lineshapes  were  normalized  so  that  g(o)  =  1  and  in 
8  (a)  the  curves  were  drawn  with  equal  linewidths.   The 
Voigt  profile  was  calculated  numerically  on  a  Hewlett- 
Packard  9310A  calculator  using  (13)  and  (24)  in  (25) .   The 
program  used  was  a  slightly  modified  version  of  the  one 
found  in  Reference  13. 
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4 .   Laser  Modes  and  Spacing 

The  general  expression  for  a  laser  resonator  cavity 
resonant  frequency  is  given  by 


V'qlm  =   q  ^  (1  +  *  +  m)  f  -21^-         (27) 


where 


,  =  ^  cos-  [(X  -  ^^,  (1  -  ^,] 


h 


R  and  R   are  the  radii  of  curvature  of  the  two  laser  mirrors 
separated  by  a  distance  L,  rj    is  the  index  of  refraction  of 
the  medium,  SL    and  m  are  the  transverse  mode  numbers  with 
values  of  0,  1,  2,.,.,    and  q  is  the  longitudinal  mode  number 
which  is  equal  to  the  number  of  half  wavelengths  in  the 
distance  L. 

From  (27)  the  separation  between  laser  oscillation 
frequencies  in  the  longitudinal  (or  Axial)  mode  for  Aq=l  is 

A^^  q  =  V(q  +  1)  -  Vq  =  V2>)L  (23) 

If  the  gain  linewidth  of  the  medium  is  greater  than  or 
equal  to  the  axial  mode  frequency  spacing  ^'^q,  several 
longitudinal  modes  may  oscillate  at  a  time. 

For  each  value  of  q  but  different  values  of  2.  +  m 
there  is  a  set  of  transverse  mode  frequencies  falling  between 
those  which  result  from  a  change  of  the  longitudinal  mode 
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index  q.   The  frequency  separation  between  adjacent  trans- 
verse modes  is 


AV  =V  -    V  =    f  — - —  (29) 

i+m         Jl+m+1  X+m  2  7^  l 


In  most  applications  of  gas  lasers,  operation  in  the  lowest 
order  ( X  =  m  =  0)  or  fundamental  transverse  mode  is  desired. 
This  can  be  accomplished  by  inserting  into  the  laser  cavity 
a  circular  aperture  of  size  such  that  it  passes  the  funda- 
mental mode  with  little  loss  but  blocks  much  of  higher  order 
transverse  modes  which  extend  further  from  the  optic  axis. 
Typical  gas  lasers  tend  to  operate  in  the  fundamental 
transverse  mode  because  of  the  aperture  represented  by  the 
discharge  tube. 

For  the  case  of  a  resonator  with  mirror  reflectivities 
^     and  ^^,  an  average  distributed  loss  constant  o(  and  the 
average  loss  per  pass  equal  to  oc  l  -  Q/y\,   (f-j^  ^  ^)      the  passive 
resonator  linewidth  is  given  by  jjRef.  14J 


^V   =   c[0C  -  4-  ^(Pl  ?  y]      ^        c  (1  -  R)     ^3Q^ 

c         2Tr'yj  2'irri   l 


where  in  the  approximation  R  =  (p  P  )   =1   and  CK  =  0 . 

J.     ^ 

In  general  the  passive  cavity  resonance  frequency  V 

does  not  coincide  with  the  molecular  line  center  frequency 

^?  ,  resulting  in  the  oscillation  at  a  frequency  near  y 

but  shifted  slightly  toward  \}    .      This  effect  is  called 

o 

frequency  pulling  and  is  due  to  a  shift  in  the  index  of 
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refraction  near  molecular  resonance.   The  atomic  dispersion 

"pulls"  the  laser  oscillation  frequency  V  from  the  passive 

resonator  value  toward  that  of  the  atomic  resonance  v) 

o 

The  frequency  of  oscillation  for  a  high  Q  resonator  cavity 
is  given  by 


V  =  V  -  (  V  -  V  )    ^^^  (31) 


where  AV  is  the  total  molecular  linewidth  due  to  Doppler 
and  pressure  broadening. 

For  operation  in  the  lowest  order  transverse  mode 
where  j2.  =  m  =  0  (2  7)  becomes 


v'q  =   [q  +  f]   ^^  =  v^  (32) 


Substitution  of  (30)  and  (32)  into  (31)  gives  for  the 
frequency  of  oscillation 


^^  J2->^L     L    47rT^2    ^2  ITrj    I,     _ 


1 
AV 


Equation  (33)  shows  that  the  laser  oscillation 
frequency  may  be  tuned  by  changing  the  cavity  length.   This 
is  the  fundamental  means  of  stabilization  techniques  using 
a  piezoelectric  driven  mirror  in  the  cavity. 


(33) 
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5 .   Frequency  Instabilities  and  Causes 

The  condition  for  stable  resonator  operation  is 
given  by  the  condition 


0  1  (1  -  Vr^)  (1  -  ^/R^)    ^     1         (34) 


The  change  in  laser  frequency  due  to  the  change  in 
cavity  length  is  given  by  [kef .  15J 


dL 
dV  =  -  V  (35) 


For  a  system  where  the  cavity  mirrors  are  fixed  to  a  cavity 

structure  with  a  coefficient  of  expansion,  cc    ,  the 

ef  f 

change  in  cavity  temperature  ^T  can  change  the  length  of 
the  cavity  resonator  and  thus  the  laser  frequency  by 


^^   -         =  -  0<  _  ^T  (36) 


V        L         eff 

Therefore  changes  in  ambient  temperature  and  acoustic 
vibrations  in  the  cavity  which  modulate  the  cavity  length 
are  sources  of  instability  which  can  contribute  to  both 
long  and  short  term  drifts  in  the  cavity  resonance  frequency. 

The  oscillation  frequency  of  a  CO   laser  can  also 
change  by  varying  functions  in  the  operating  parameters 
such  as  pressure  and  discharge  current.   Mocker  in  Reference  16 
reported  a  5  -  3  MHz/torr  downshift  in  frequency  due  to  an 
increase  in  pressure  and  a  .5  -  .9  MHz/torr  upshift  in 
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frequency  due  to  an  increase  in  discharge  current.   He 
concluded  that  the  pressure  dependent  frequency  shift  is 
caused  by  a  change  in  the  overall  index  of  refraction  of 
the  gas  discharge  and  the  current  dependent  frequency 
shift  is  attributed  to  the  change  in  the  index  of  refraction 
caused  by  the  change  in  the  electron  density  of  the  plasma. 
In  a  CO  gas  discharge  the  refractive  index  -rj   is  determined 
by  all  constituents  in  the  plasma  such  as  the  neutral 
molecules  and  electrons:  [kef.  17j 

2 
"^   =    1   -    h   — --^ T    +  27rCK     r?         (37) 

Where  V      and  M  are  the  electron  density  and  mass 
'e      e 

respectively:  Yj        and  oc   are  the  density  and  polarizability 
of  the  neutral  molecules  respectively.   A  significant  re- 
duction of  the  refractive  index  may  arise  from  dissociation 
of  CO2  molecules  into  components  with  a  smaller  electric 
permittivity. 

Therefore  it  is  important  to  operate  a  CO   laser 
with  a  well  regulated  low  ripple  power  supply  and  closely 
controlled  gas  pressure  and  temperature. 

The  frequency  deviation  ^V  from  the  resonant 
frequency  V  is  therefore  determined  by  the  fluctuation  in 
the  refractive  index  6^-f^   and  cavity  length  ^L  and  is  given  by 

^^   =   ^^   +  -^22-  (38) 
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For  a  sealed  CO   laser  the  pressure  fluctuations  are 
negligible  and  the  temperature  and  flow  rate  of  cooling 
water  are  the  dominant  causes  of  instability  in  the 
laser . 
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B.   FREQUENCY  TUNING  MECHANISM 

By  bringing  a  cavity  mode  into  resonance  with  a 

vibrational-rotational  peak  a  CO   laser  can  be  tuned  to  a 

vibrational-rotational  line.   The  cavity  mode  spacing  (fv 

of  a  single  mode  TEM    laser  dictates  the  chance  of 

ooq 

oscillation  on  a  specific  rotational  level  and  changes 
as  the  cavity  mirror  separation  is  varied.  The  optical 
length  of  the  Model  941  laser  cavity  is  45  cm,  short 
enough  to  ensure  single  axial-mode  operation  and  give  a 
value  for  6^c=  "^q  -  "^q+l  =  V2L  of  333MHz.  The  cavity 
linewidth  calculated  from  (30)  with  c<  =  o  and  using  the 
reflectivities  given  by  Sylvania  for  the  cavity  mirrors, 

f  =  1.0  and  p2  =  -"^^  ^^  19MHz.   Each  vibrational- 
rotational  transition  is  lineshape  broadened  with  a  line- 
width  AV  .   From  (2  3)  the  Doppler-broadened  linewidth 
for  the  CO   transitions  is  approximately  50MHz.   For  the 
pressures  and  mixtures  used  in  single  transverse  mode 
(small  bore  diameter)  CO   lasers  (15  -  20  Torr)  the 
Lorentzian  pressure  broadening  ^v',_  can  exceed  the  Doppler 
broadening.   The  value  listed  by  Sylvania  for  ^Vj_  is  lOOMHz 
giving  a  total  linewidth  using  (26)  of  112MHz.   The  absolute 
frequency  separation  of  the  vibrational-rotational  lines 

AV     is  between  52  and  57GHz  depending  on  the  J  number 
v-R 

of  the  transition.   Figure  9,  although  not  drawn  to  scale, 
shows  three  P(J)  lines  and  a  set  of  cavity  modes  showing 
the  frequency  spacing  involved.   The  successive  order  of 
the  transitions  is  dictated  by  the  absolute  cavity  length; 
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they  occur  in  an  arbitrary  fashion  giving  each  laser  its 
own  spectral  signature. 

As  a  CO   laser  is  tuned  in  frequency  by  changing  the 
cavity  length,  competition  effects  occur  and  oscillation 
switches  from  one  rotational  line  to  another.   In  the 
interrotational  level  competition  region,  two  rotational 
levels  compete  for  the  total  population  of  the  same  and 
higher  vibrational  levels.   Mocker  [^Ref.  is]  reported  that 
the  competition  effects  are  due  to  intensity-dependent 
anomalous  dispersion  arising  from  index  and  gain  saturation 
and  that  the  switching  from  one  rotational  level  to  an 
adjacent  level  occurs  over  a  small  finite  frequency  interval. 

While  the  laser  oscillates  on  a  given  vibrational- 
rotational  transition  the  population  density  or  the  popu- 
lation inversion  on  that  particular  transition  is  reduced 
by  saturation  effects  until  the  single  pass  gain  in  the 
medium  equals  the  single  pass  loss  of  the  optical  cavity. 
If  the  vibrational  level  lifetime  (typical  values  1  -  10msec) 
is  considerably  longer  than  the  rotational  relaxation  time 
{<  1  jm^sec)    a  Boltzmann  distribution  will  describe  the 
population  density  of  the  various  rotational  levels  at 
the  same  temperature.   Therefore  if  the  gain  on  the 
oscillating  transition  is  the  highest  the  gain  on  the 
non-oscillating  vibrational-rotational  transitions  starting 
from  other  rotational  levels  in  the  vibrational  level  will 
be  less  than  the  single-pass  loss  of  the  optical  cavity. 
This  means  that  once  the  laser  transition  starts 
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oscillating,  the  gain  on  the  other  transitions  falls  below 
that  required  for  laser  oscillation  with  the  other  rotational 
levels  coupling  (due  to  collisions)  their  stored  energy 
into  the  lasing  transition  and  only  one  transition 
oscillates.   This  is  the  situation  of  very  strong  com- 
petition effects.   As  the  cavity  mode  is  tuned  away  from 
the  center  of  the  lasing  transition  a  point  is  reached 
where  another  transition  has  equal  gain  and  a  competition 
will  be  initiated.   The  intensity  dependent  dispersion 
represents  a  frequency  pulling  term  such  that  a  cavity 
resonance  is  pulled  toward  the  line  center  frequency. 
Transitions  of  approximately  equal  strength  have  their 
competition  region  halfway  between  line  centers,  while 
unequal  strength  transitions  compete  closer  toward  the 
weaker  transition.   Figure  10  illustrates  the  competition 
between  vibrational-rotational  transitions  and  how  the 
laser  can  be  tuned  in  frequency  by  varying  the  cavity  mode 
separation  as  a  result  of  changing  the  cavity  length.   In 
Figure  10(a)  the  cavity  mode  at  V^^' ,  is  in  resonance  with 

the  P(J,)  laser  transition  causing  it  to  oscillate  and  P(J^) 
1  ^ 

is  between  the  cavity  modes.   In  Figure  10(b)  the  cavity 
length  is  increased  slightly  and  the  modes  at  \)^^   and  y^g. 
are  in  resonance  with  the  P(J^)  and  P ( J  )  transitions 
respectively  and  competition  can  occur.   If  the  gains  are 
mutually  exclusive  then  only  one  line  would  survive  at  the 
expense  of  the  other.   If  the  gain  of  each  line  is  sufficient 
to  sustain  oscillation  simultaneously  then  both  laser  lines 
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would  coexist.   In  Figure  10(c)  the  cavity  mode  at  V^,  is 
in  resonance  with  the  P(J-)  transition  causing  it  to  oscillate 
while  the  P ( J  )  transition  is  between  the  cavity  modes.   This 
frequency  tuning  mechanism  makes  it  possible  to  use  a  piezo- 
electric driver  to  actively  control  and  stabilize  the 
frequency  output  of  the  laser. 
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3  CO-  LASER  P(J)  TEANSITIONS  AND  CAVITY  MODES 
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C.   FREQUENCY  STABILIZATION  TECHNIQUES 

Although  many  of  the  instabilities  contributing  to 
frequency  fluctuations  in  the  CO   laser  can  be  reduced  by 
external  means  (ie.  control  of  discharge  current  by 
regulated  power  supply  and  control  of  system  temperature 
with  the  aid  of  a  separate  temperature  controlled  heat 
exchanger)  to  frequency  stabilize  the  laser  actively 
requires  the  ability  to  change  either  the  index  of  refraction 
of  the  laser  medium  or  resonator  cavity  length.   The  laser 
cavity  length  can  be  accurately  controlled  by  mounting  a 
piezoelectric  crystal  which  expands  and  contracts  proportion- 
ally to  an  applied  voltage.   An  error  or  frequency  correction 
signal  that  is  proportional  to  and  indicates  the  direction 
of  the  required  length  correction  can  drive  the  piezoelectric 
crystal  and  frequency  stabilize  the  laser.   A  common  feature 
shared  by  a  number  of  approaches  involves  the  basic  principle 
of  establishing  a  frequency  reference  either  internal  or 
external  to  the  laser  to  which  the  frequency  is  locked. 
By  applying  a  modulation  signal  to  the  laser,  causing  its 
oscillating  frequency  to  "dither'  about  a  reference  point, 
and  using  phase  sensitive  detection  techniques,  an  error 
signal  can  be  fed  back  to  servo  the  laser  frequency  to  the 
reference. 

Another  stabilization  technique  which  has  some  useful 
advantages  over  the  methods  referred  to  above  uses  a  frequency 
reference  that  is  external  to  the  laser  with  the  discriminant 
obtained  not  by  dither  modulating  the  laser  frequency  about 
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the  reference,  but  by  modulating  the  reference  about  the 
laser  frequency.   The  important  advantage  of  this  technique 
is  that  it  avoids  impressing  modulation  (AM  or  FM)  on  the 
laser  output.   A  gas  cell  containing  a  molecular  absorber 
which  is  tuned  into  resonance  with  the  desired  CO   laser 
line  by  a  suitable  D.C.  electric  field  applied  to  the 
absorber  provides  the  frequency  reference.   By  super- 
position of  an  A.C.  electric  field,  the  absorption  line 
center  can  be  modulated  about  the  lasing  frequency  and 
phase  detected,  producing  the  required  correction  signal 
back  to  the  piezoelectric  drive  on  the  laser.   Another 
advantage  of  this  technique  is  the  fact  that  it  is  voltage- 
controllable  which  means  the  laser  frequency  may  be  con- 
trollably  tuned  across  the  transition  line  width  by 
programming  the  D.C.  bias  voltage  applied  to  the  Stark 
cell. 
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D.   STARK  CELL  STABILIZATION 
1.   Stark  Effect 

Stark  effects  are  the  changes  in  the  spectrum  of 
a  system  due  to  an  applied  electric  field.   The  rotational 
spectrum  of  a  molecule  which  has  a  dipole  moment  is  modified 
in  an  electric  field,  since  the  field  exerts  torques  on  the 
electric  dipole  moment  changing  its  rotational  motion. 

The  effect  of  an  electric  field  on  the  rotational 
energy  levels  of  molecules  is  given  by  quantum-mechanical 
perturbation  theory.   The  Stark  effect  in  a  gas  occurs  when 
an  external  electric  field  is  applied  to  the  gas.   The  effect 
is  characterized  by  the  splitting  of  an  absorption  peak 
into  several  peaks  by  removing  the  degeneracy  of  the 
rotational  energy  levels.   The  application  of  the  field 
in  the  absorbing  gas  results  in  the  splitting  of  a  line 
into  2J  +  1  lines . 

The  first  order  Stark  effect  which  is  linear  with 
the  applied  electric  field  is  produced  by  applying  the  elec- 
tric field  to  a  symmetric  top  molecule  (molecule  with  three 
or  more  fold  symmetry  -  molecule  with  two  principal  moments 
of  inertia  equal).   A  symmetric  top  molecule's  principal 
advantage  is  the  existence  of  a  permanent  electric  dipole 
moment.   The  Stark  splitting  of  the  energy  levels  in  molecules 
possessing  a  permanent  dipole  moment  is  given  by  [Ref.  19] 


A(  M  E  K  ,   ^ 

hAV  =  -^^-- r-   cos  (9-  (39 

J  J  +  1) 
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Where  >*-(  is  the  molecular  dipole  moment  and  E  is  the  electric 
field  strength,  <3  is  the  angle  between  >\  and  E,  J  is  the 
total  rotational  quantum  number,  M  is  the  magnetic  quantum 
number  or  the  component  of  J  along  a  fixed  direction,  and 
K  is  the  quantum  number  describing  the  interaction  of  the 
two  principal  moments  of  inertia.   The  selection  rules  for 
vibration-rotational  transitions  are 

AJ  =  0,  ±1 

AK  =  0   For  Parallel  Bands 

AK  =  +1   For  Perpendicular  Bands 

AM  =  0   For  Stark  Field  Parallel  to  Laser 
Electric  Field 

AM  =  ±1  For  Stark  Field  Perpendicular  to 
Laser  Electric  Field 

For  asymmetric  molecules  or  those  which  do  not 
have  a  permanent  electric  dipole  moment,  the  expression 
for  the  Stark  splitting  of  the  energy  levels  is  more 
complicated  and  involves  second  order  perturbation  calcu- 
lations  including  the  neighboring  states  J   =  J  +  1  and 
J  =  J  -  1.   The  energy  splitting  goes  as  the  square  of  the 
external  electric  field.   The  second  order  quadratic  Stark 
effect  is  usually  considered  negligible  in  Stark  modulation. 
[Ref.  20] 

For  slightly  asymmetric  top  molecules  an  inter- 
mediate type  of  Stark  effect  results.   These  molecules 
can  have  pairs  of  nearly  degenerate  levels  and  a  first 
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order  Stark  effect  can  result  for  small  molecular 
asymmetry  and  large  fields. 
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2 .   Selection  of  Gases 

Changes  in  absorption  produced  by  an  applied 
electric  field  on  the  states  involved  in  transitions  as 
given  by  the  vibrational-rotational  selection  rules  may 
result  in  a  modulation  of  the  laser  intensity  by  a  gas 
containing  an  appropriate  molecule.   In  order  to  achieve 
an  appreciable  Stark  effect  and  large  modulation  J  should 
be  small. 

The  selection  of  a  suitable  molecule  for  modulation 
of  an  individual  laser  is  rather  straightforward.   The 
criteria  to  achieve  a  high  modulation  index  are:   1)  it  is 
necessary  to  use  molecules  with  a  large  dipole  moment; 
2)  slightly  asymmetric  molecules  may  be  used  provided  the 
Stark  effect  is  essentially  linear  at  the  values  of  applied 
electric  field  used  for  modulation;   3)  a  linear  Stark 
effect  in  either  lower  or  upper  state  or  both  is  suitable 
and  the  overall  effect  is  determined  by  combination  of  the 
Stark  effect  on  both  states;   4)  the  absorbing  Stark  sensi- 
tive line  and  laser  line  should  be  within  one  or  two 
Doppler  widths  of  the  absorbing  gas.   This  requirement  is 
not  valid  in  the  pressure  broadening  case  where  the  higher 
pressures  result  in  broader  lines,  poorer  resolution  and  a 
smaller  index  of  modulation  with  equivalent  fields. 

Investigations  by  Landman,  Marante  and  Early  JRef.  2lJ , 
Jensen  and  Tobin  [Ref.  2  2j  and  most  recently  Martin,  Corcoran 
and  Smith  [Ref.  2  3]  have  reported  many  gases  capable  of 
Stark  modulation  of  a  CO   laser.   Table  II  lists  some  of 
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the  absorbing  gases  used  for  Stark  modulation  and  their 
corresponding  CO   laser  lines. 


TABLE  II 

GASES  USED  IN  STARK  MODULATION  OF  00° 1  -  10°0  TRANSITION 
Gas  Laser  Line  Modulated 


C2H^C1  P(16),  P(18),  P{20),  R(16),  R(18),  R(20) 

CH^Br  P(22) 

CH^I  P(14) 

NH  R(14) 

C2H4F2  P(14),  P(18),  P(20),  R(18) 

C2H2F2  P(10),  P(12),  P(14),  P(18),  P(20),  P(22) 

NH2D  P(14) ,  P(20) ,  R(12) 
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3 .   NH2D  Rotational-Vibrational  Spectra 

Of  the  gases  in  Table  II  the  microwave  and  infrared 
spectra  of  NH^D  have  been  studied  extensively.   Brewer  et .  al. 
in  Reference  24  reported  a  rapid  Stark  tuning  rate  and 
relatively  high  absorption  coefficient  for  NH  D  making  it 
possible  to  obtain  substantial  modulation  with  relatively 
low  electric  field  strength  and  short  path  length.   Kelly, 
Francke  and  Feld  [Ref.  25j  observed  the  transitions  which 
could  be  tuned  into  resonance  with  three  of  the  10 . 6^m  CO2 
laser  lines,  namely  P(20),  P(14)  and  R(12).   The  NH  D 
transitions  were  identified  as  belonging  to  the  type 

(v-,  J)  =  (0,  4)* *(1,  5).   The  assigned  transitions  are 

shown  in  Figure  11  with  the  corresponding  CO   laser  lines 
enclosed  in  brackets.   The  lower  and  upper  levels,  the  Stark 
shift  (S      required  to  bring  each  NH  D  transition  into  resonance 
with  its  respective  laser  line,  laser  line  wavelength  (cm~-^) 
and  zero  field  splitting  (cm~  )  for  the  three  laser  lines 
are  given  in  Table  III  from  Reference  24.   Of  primary 
interest  in  this  discussion  is  the  transition  involved  in 
absorption  of  the  P(20)  CO^  laser  line  in  the  10.4/^m  band. 

The  NH^D  absorption  line  involved  belongs  to  the 
transition 


^^'  ^k/k,  ^  =  ^^a'  ^04)  ^--*   (la'  ^05^ 


This  transition  is  nearly  in  coincidence  with  but  at  a 
slightly  lower  frequency  than  the  944.19cm~   P(20)  line 
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of  the  CO2  laser.   The  P(20)  laser  frequency  exceeds  the 
zero  field  splitting  of  the  (O  ,  4   )  «* — m{l    ,    5-^  ) 
transition  by  1720  ±  60  MHz.   A  relatively  large  and  linear 
Stark  effect  results  from  mixing  of  the  two  nearly  degenerate 

rotation-inversion  levels  (O  ,  ^csa)    ^^*^  ^^s'  ^^A^    ^^   *'^® 
lower  state  of  the  transition.   The  subscript  s(a)  refers 
to  a  symmetric  (antisymmetric)  inversion  state  with  respect 
to  simultaneous  exchange  of  all  pairs  of  identical  nuclei. 
The  upper  level  (1  ,  5q_)  has  only  a  negligible  second 


order  Stark  effect  [Ref .  26j  . 
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RELEVANT  ElfiilRGY  LEVELS  OF  NH2D 
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TABLE  III 

IDENTIFICATION 

OF  OBSERVED  NH^D  TRANSITIONS 

LASER 
LINE 

LOWER 
LEVEL 

UPPER 
LEVEL 

(MHz) 

LASER  FREQ, 
(cm-') 

ZERO  FIELD 
SPLITTING  (cm^') 

P(14) 

°a»  V 

^s'  hu 

-2300±60 

949.43 

949.40 

P(20) 

°a'  V 

\>   ^05 

-1720*60 

94^.18 

944.14 

R(12) 

0.»  ^1/ 

T  »  5,/ 

+2800i1 00 

970.55 

970.64 
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4.   Error  Signal  Analysis 

The  modulated  output  intensity  from  a  Stark  cell 
due  to  absorption  by  a  single  Stark  component,  exhibiting 
a  linear  Stark  effect  is  given  by 


I  =  Iq  exp  -['^i]  (40) 


Here  T   is  the  absorption  coefficient  of  the  gas  as  a  function 

of  light  intensity  and  applied  electric  field,  I   is  the 

o 

input  intensity  and  £  is  the  cell  length.   The  line  center 
absorption  coefficient  where  the  line  is  collision-broadened 
in  the  high  pressure  Lorentzian  region  is  given  by  JRef .  27j 

^ '  ^»'  =  -  ^  ^  '  T^  '  ("z  -  It  V  (41) 

Where  Xq is  its  wavelength,  ^^     is  the  radiative  lifetime 

of  the  vibrational-rotational  transition,  N   and  N   are  the 

upper-  and  lower-level  population,  respectively,  g^   and  g 

^  1 

are  the  corresponding  statistical  weights  and  ^V  is  the 

homogeneous  linewidth  as  given  in  (19) .   The  pressure 

dependence  of  ^V,  is  cancelled  since  N   and  N   are  also 

^  2      1 

proportional  to  pressure,  and  A(?v.o)  is  a  constant. 

The  absorption  coefficient  at  line  center  for  an 
inhomogeneous  line  in  the  low-pressure  Doppler  broadened 
limit  is  given  by  (Ref.  2  7j 

MXc,)    = -rr^ i-^ — )   (N^ ^  N  )     42) 

8Tr  ^21    ^Vjj      ^        2    gj_    1 
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Where  AV^^  is  the  inhomogeneous  linewidth  given  by  (2  3)  and 

is  a  constant.   The  pressure  dependence  of  N   and  N   cause 

2       1 

A (  Xo  )  to  vary  linearly  with  pressure. 

Plant  and  Abrams  in  Reference  27  plotted  the  line 
center  absorption  coefficient  vs  Stark  cell  pressure  for 
NH  D  as  shown  in  Figure  12.   Using  Voigt  profile  fits  to 
the  line  shapes  they  were  able  to  get  a  relation  for  the 
line  center  absorption  coefficient  which  was  an  excellent 
fit  to  the  experimental  points.   This  relation  is 

'^^  ^°'  =  -  W  =Ffl  -Jj^^    'V  if-"!'  '1--=^^^'  (expCa^3)     (43) 

Where  a  =  (  ^"^l/ ^V^^  )  {  JbY\2)       and  erf  is  the  error  function. 
From  Figure  12  at  pressures  above  about  5  torr  one 

can  approximate  the  Stark  cell  absorption  lineshape  with  a 

-1 
Lorentzian  and  A (  X^  )  =  0.023cm   ,  a  constant.   The  mod- 
ulator bandwidth  reaches  a  limit,  since  the  Stark  effect 
changes  character  as  the  frequency  of  the  perturbing  field 
becomes  comparable  to  the  linewidth.   Claspy  et .  al.  [Ref.  28J 
estimated  that  the  modulator  bandwidth  increased  with 
pressure  according  to  a  pressure  broadening  coefficient. 
Therefore  the  broader  line  at  higher  pressures  required  a 
greater  field  variation  for  the  same  absorption  change. 

The  absorption  coefficient  in  (40)  can  be  given  by 


r  =  A(  K)^^^   «  g(y  )  (44) 

MAX 


53 


LINE  -  CENTER  ABSORPTION  COEFFICIENT  VS  STARK  CELL  PRESSURE 
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Where  oc  is  a  constant  dependent  on  gas  pressure  and  g(V  ) 
is  the  lineshape  function  of  the  absorbing  gas.   Assuming 
a  Lorentzian  lineshape  and  letting  c<  equal  the  pressure 
broadening  coefficient  C,    times  the  pressure  P,  (44)  can 
be  written  as 


r=A(Xo)     CP  = ^^n    ^^ -5-=-  (45) 

MAX      2Tr[(V-V,)^-(^V/2y] 


where  A "^^  is  the  linewidth  of  the  Stark  cell,  and  (  V-Vq  ) 
is  the  frequency  offset  of  the  absorption  line  center  from 
the  laser  frequency.   Let  (  V-Vq)  be  defined  as  6  where 


€  =  €.   +  SE  (46) 

o 


€   is  the  offset  frequency  dependent  on  the  offset  of  the 
absorption  line  from  the  laser  frequency  at  zero  field  and 
S  is  the  Stark  coefficient  or  Stark  tuning  rate  which  relates 
the  line  shift  to  the  applied  electric  field. 
Substituting  (46)  into  (45)  gives 

r  =  A(  Xo  )     CP  = ^^ ^r-   (47) 

MAX      27r[(€  ^  +  SE)^  +  (  ^V/2)^J 

Therefore,  for  a  fixed  pressure  in  the  cell  and  a  set  laser 
line,  T   the  absorption  coefficient  is  only  dependent  on  the 
applied  Field  E  .   Impression  of  an  A.C.  voltage  with  an 
electric  field  in  the  form  E  =  AE  cosojt  on  the  applied 
field  E  will  give  (T  a  time  dependence  by  oscillating  the 
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absorption  line  shape  about  the  line  center  frequency 
which  will  act  as  a  dither  mechanism  about  V  .   A  schemiatic 
of  the  dither  mechanism  for  a  Stark  cell  modulator  is  shown 
in  Figure  13.   The  remainder  of  the  work  described  in  this 
thesis  will  refer  to  the  P(20)  -  NH  D  combination. 

Nussmeier  and  Abrams  [Ref.  263  give  the  frequency 
offset  between  the  absorption  line  center  of  NH  D  and  the 
P(20)  10.6  jL^m   line  center  as  a  function  of  applied  voltage  as 


€  =  e   +  SE  =  -2042  +  0.143|M|E  (48) 

o 


V^ere  €.  is  given  in  MHz,  M  is  the  e  component  of  the 

rotational  quantum  number  J  and  has  integer  values  from 

-4  to  4,  and  E  is  the  applied  electric  field  in  /cm. 

Setting   6=0  and  M  =  ±4  in  (48)  shows  that  an  electric 

field  of  3570   /cm  is  needed  to  bring  one  component  of 

NH  D  absorption  into  resonance  with  the  P(20)  transition. 
2 

The  values  reported  in  the  literature  for  the  Stark  tuning 
rate,  S,  for  the  P(20)  CO   line  and  NH  D  with  M  =  ±4  are 
consistent  with  Nussmeier  and  Abrams  value  of  .572  MHz   /v 
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5 .   Stark  Cell  Stabilization  System 

A  schematic  of  a  laser  stabilization  system  is 
shown  in  Figure  14  and  is  a  somewhat  simplified  version  of 

the  one  given  in  Reference  29.   A  small  percentage  (less 

Iw ,   2 
than   /cm  to  avoid  saturation  effects)  of  the  laser 

output  is  directed  into  the  Stark  cell  via  a  beam  splitter 

and  plane  mirror.   The  basic  Stark  cell  consists  of  two 

parallel  metal  electrodes  in  a  vacuum  tight  container  with 

an  entrance  and  exit  window.   A  precision  high-voltage 

power  supply  biases  one  Stark  electrode  to  bring  an 

absorption  line  in  NH  D  into  resonance  with  the  P(20) 

line  center.   An  A.C.  voltage  from  an  audio  oscillator 

drives  the  second  electrode  at  a  set  frequency.   The  cell 

contains  NH  D  formed  by  mixing  one  part  NH  with  ND   to  a 

desired  pressure  in  the  cell.   The  exchange  of  hydrogen 

and  deuterium  is  reported  in  Reference  30  to  be  rapid  with 

an  equilibrium  concentration  to  45%  NH  D  reached  in  less 

than  one  minute.   When  the  laser  is  oscillating  on  the 

P(20)  line  the  Stark  effect  will  produce  a  sinusoidal 

modulation  in  the  intensity  through  the  cell  which  can  be 

detected  by  a  pyroelectric  or  HgCdTe  detector.   The  signal 

is  phase  sensitively  detected  relative  to  the  reference 

audio  oscillator  waveform  in  a  lock-in  amplifier.   The 

detected  phase  is  dependent  on  the  relative  frequencies  of 

the  absorption  line  center  and  the  laser  output.   The 

magnitude  and  sign  of  the  lock-in  amplifier  D.G.  output 

which  is  proportional  to  the  phase  difference  between  the 
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laser  frequency  and  the  Stark-shifted  reference  will  produce 
a  suitable  error  signal.   After  this  error  signal  is 
appropriately  phased,  integrated  and  amplified  it  is  fed 
back  to  the  piezoelectric  mounted  laser  cavity  mirror  to 
tune  out  the  frequency  difference  and  the  frequency  control 
feedback  system  is  complete.   Because  the  laser  is  capable 
of  oscillating  on  a  number  of  vibrational-rotational  laser 
transitions  dependent  on  the  resonator  cavity  length,  a 
linear  ramp  voltage  applied  to  the  piezoelectric  drive  will 
cause  the  laser  to  scan  through  its  spectral  signature 
profile  until  the  P(20)  transition  is  reached  at  wnicn  time 
tne  servo  stabilization  technique  can  be  initiated. 
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III.   EXPERIMENTAL  INVESTIGATION 

Preliminary  engineering  work  done  in  Reference  13 

studied  the  characteristics  of  the  Sylvania  CO  Model  941 

2 

Laser  pertinent  to  the  design  of  a  Stark  cell  stabilized 
system.   The  instabilities  of  the  laser  due  to  fluctuations 
in  discharge  current  caused  by  an  unregulated,  unfiltered 
laser  power  supply  and  system  temperature  variations  were 
measured  and  discussed.   The  piezoelectric  parameter  and 
thermal  expansion  coefficient  also  were  determined  for  the 
laser.   Computer  analysis  demonstrated  the  feasibility  of 
producing  a  suitable  error  signal  for  phase  sensitive 
detection  with  available  equipment  and  the  need  for  a  servo 
loop  integrator.   A  vacuum  tight  Stark  cell  was  constructed 
and  tested  to  3000  volts  with  no  evidence  of  breakdown. 

The  parameters  for  the  Sylvania  941  piezoelectric 
tunable  and  linearly  polarized  CO-  laser  used  in  the 
following  discussion  of  experimental  procedures  and  results 
are  given  in  Appendix  A.   Appendix  B  discusses  the  optical 
design  and  specifications  of  the  optical  engineering  CO 
laser  spectrum  analyzer  used  in  the  experiments. 

To  prevent  the  instabilities  caused  by  changes  in  dis- 
charge current  a  stable  high  voltage  power  supply  (Sylvania 
941S)  which  in  turn  was  run  from  a  stable  line  regulator 
was  used  to  power  the  laser.   The  use  of  a  circulation  pump 
minimized  flow  rate  fluctuations  and  the  system  usually 
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reached  an  equilibrium  temperature  condition  after  a  few 
minutes.   The  current  to  the  laser  from  the  941S  power 
supply  could  not  be  controlled  manually  and  was  preset 
for  optimum  performance  of  the  laser. 
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A.   SERVO  LOOP  INTEGRATOR 

In  order  to  provide  optimum  frequency  stabilization 
of  the  CO   laser  the  P(20)  line  should  be  maintained  in 
coincidence  with  the  center  absorption  frequency  of  a 
Stark  shifted  NH  D  transition  as  accurately  as  possible. 
The  two  transitions  in  coincidence  will  produce  a  zero 
error  signal.   As  the  P(20)  line  drifts  in  frequency  the 
magnitude  and  sign  of  error  signal  will  be  proportional  to 
the  change  in  frequency,  and  direction  of  the  drift, 
respectively.   If  a  proportional  correction  signal  is 
supplied  directly  to  the  piezoelectric  driver  the  magnitude 
of  the  error  signal  will  decrease  to  an  equilibrium  condition 
at  a  nonzero  value  for  € .      The  equilibrium  condition  exists 
when  the  P3  voltage  producing  a  frequency  difference  e  is 
equal  to  the  P3  voltage  supplied  by  the  correction  signal. 

Integration  of  the  error  signal  will  produce  a  correction 
signal  that  will  drive  the  piezoelectric  crystal  in  a  manner 
proportional  to  the  total  time  integral  of  the  error  signal. 
This  integration  will  allow  the  error  signal  to  return  to 
zero  at  the  required  P2  voltage.   The  integrated  signal 
will  continue  to  increase  until  the  frequency  difference 
(error  signal)  is  zero.   At  £  =  0  the  integrator  maintains 
the  voltage  at  a  constant  D.C.  value  until  a  nonzero  error 
signal  indicates  the  need  for  a  change. 

The  integrator  shown  in  Figure  14  serves  a  dual  purpose: 
i)  integration  of  the  error  signal  supplied  by  phase  detection 
in  the  lock-in  amplifier  when  the  stabilization  feedback 


63 


control  system  is  in  operation.   ii)  Integration  of  a 
selectable  D.C.  voltage  to  provide  a  linear  ramp  voltage 
to  the  laser  piezoelectric  driver  in  order  to  search  for 
the  required  laser  transition  before  the  servo  stabilization 
procedure  can  be  initiated. 

The  integrator  described  in  Appendix  C  was  designed 
to  respond  to  both  slow  and  rapid  changes  in  frequency  of 
the  P(20)  line  from  the  Stark  shifted  NH  D  absorption  center 
frequency.   The  gain  was  designed  to  be  equal  to  one  to  make 
the  integrator  output  equal  to  the  integral  of  the  input. 

A  diagram  of  the  apparatus  used  to  examine  the 
characteristics  of  the  integrator  is  shown  in  Figure  15 . 

The  input  to  the  integrator  was  provided  by  a  regulated 
0  -  25  V  D.C.  supply  which  was  set  at  a  constant  voltage 
A  three  way  reversing  polarity  switch  S:  was  used  to  provide 
positive  or  negative  values  of  V;_  to  the  integrator  or  to 
isolate  the  integrator  from  the  input  voltage.   The  input 
voltage  to  the  integrator  was  fed  into  channel  A  of  a 
time  base  plotter  and  the  output  of  the  integrator  into 
channel  B.   Figure  15  shows  the  integrator  output  as  a 
function  of  a  step  input  with  V|^  =  ±2    volts.   The  graphical 
output  indicates  that  integration  of  a  2  volt  input  for 
approximately  10  seconds  gives  an  integration  output  of 
20  volts  and  slope  of  1.95   /sec.   Thus  the  integrator 
output  equals  the  integral  of  the  input.   Figures  17  a  &  b 
respectively  show  typical  operation  and  response  of  the 
integrator  to  step  inputs  of  ±4V  and  ±6V.   Since  the 
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integrator  is  part  of  the  feedback  control  system  that 
drives  the  piezoelectric  crystal,  differences  between  the 
analog  integration  of  the  input  signal  and  the  exact  math- 
metical  values  are  not  critical  and  can  be  compensated  for 
in  the  proportionality  constant  associated  with  amplification 
of  the  correction  signal.   In  other  words,  the  integrator 
is  concerned  primarily  with  the  manner  in  which  the  piezo- 
electric crystal  is  driven  rather  than  the  magnitude  of  the 
correction  signal. 

The  integrator  saturates  at  a  voltage  of  +14  volts; 
therefore  if  the  P(20)  line  consistently  drifts  in  the  same 
direction  the  integrator  may  saturate.   By  ganging  a  polarity 
reversing  switch  to  the  input  and  output  of  the  integrator, 
saturation  can  be  controlled  by  reversing  the  sign  of  the 
signal  in  magnitude  while  maintaining  sign  of  the  proper 
correction  signal  to  tune  out  the  frequency  difference. 
Minor  drift  in  steady  state  voltage  output  inherent  in 
analog  integrations  can  be  controlled  by  a  voltage  offset 
null  on  the  integrator.   This  possible  integrator-induced 
long  term  instability  will  cause  the  laser  frequency  to 
drift  slightly  but  it  can  be  corrected  for  in  the  feedback 
control  system  like  any  other  instabilities  in  the  system. 

The  incorporation  of  an  integrator  should  improve  the 
ability  of  the  feedback  control  system  to  frequency  stabilize 
the  laser  and  produce  a  system  suitable  for  laboratory  or 
field  use. 
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B.   LASER  SPECTRAL  SIGNATURE 

The  laser  is  capable  of  oscillation  on  a  number  of 
transitions  other  than  the  P(20)  as  the  resonator  optical 
length  varies.   This  capability  adds  to  the  versatility  of 
the  laser  as  well  as  satisfying  a  need  to  search  for  and 
find  the  desired  transition  before  the  feedback  stabilization 
system  can  be  initiated.   If  the  cavity  length  of  a  CO2 
laser  is  slowly  increased  a  certain  sequence  of  lines  will 
appear.   After  changing  the  cavity  length  approximately  one 
half  wavelength  the  sequence  will  repeat  itself.   The  output 
profile  obtained  by  changing  the  laser  mirrors  separation  is 
called  the  laser  signature. 

A  schematic  of  the  experiment  to  determine  the  spectral 
signature  of  the  laser  is  shown  in  Figure  13.   The  spectral 
signature  was  produced  by  plotting  the  relative  intensity 
of  the  laser  as  a  function  of  applied  voltage  to  the 
piezoelectric  driver  (PZ  voltage) .   The  P3  voltage  was 
supplied  by  a  Kepco  model  ABC  regulated  D.C.  supply  with  a 
range  of  0  -  1500  volts.   It  was  operated  in  the  voltage- 
programmable  mode  with  a  positive  ground.   It  was  pro- 
grammed to  operate  at  a  voltage  ratio  of  100:1  output  to 
input.   Using  the  integrator  described  in  Appendix  C  to 
supply  a  linear  ramp  voltage  as  the  input  to  the  programmable 
power  supply,  the  laser  could  be  scanned  through  a  free 
spectral  range.   The  scan  rate  was  controlled  by  the  magni- 
tude of  the  input  voltage  to  the  integrator  since  the  slope 
of  the  integrator  output  voltage  is  directly  proportional 
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to  the  input  voltage.   The  direction  of  scan  was  controlled 
by  a  3  position  reversing  polarity  switch  that  changed  the 
sign  of  the  slope  as  a  result  of  changing  the  sign  of  the 
input  voltage.   The  output  of  the  integrator  which  was 
directly  proportional  to  the  P3  voltage  was  supplied  to  the 
X  axis  of  an  X  -  y  plotter.   The  relative  intensity  of  the 
chopped  laser  output  was  obtained  from  a  pyroelectric 
detector  signal  that  was  amplified,  filtered  and  rectified. 
The  relative  intensity  signal  was  then  supplied  to  the 
y  channel  of  the  plotter.   The  frequency  of  the  laser  for 
a  given  PS  voltage  could  not  be  monitored,  but  the 
vibrational-rotational  transition  could  be  identified 
using  an  Optical  Engineering  CO   spectrum  analyzer  described 
in  Appendix  B.   Typical  results  are  shown  in  Figure  19  for 
an  increasing  P3  voltage  scan  (0  -  1500  V  D.C.)  and  Figure  20 
for  a  decreasing  Pg  voltage  scan  (1500  -  0  V  D.C.) . 

The  power  meter  in  Figure  18  monitored  the  operation 
of  the  laser.   It  was  used  as  a  comparison  with  the  results 
shown  in  Figures  19  &  20.   Figures  21  &  22  give  typical 
results  for  relative  power  output  as  a  function  of  applied 
P2  voltage. 

In  each  set  of  plots  the  scan  direction  appears  to 
have  caused  a  hysteresis  effect  [Ref .  24]  .   Differences  in 
the  traces  appear  to  be  caused  by  hysteresis  effects  and 
nonlinearity  of  the  piezoelectric  element.   If  the  dipoles 
(whatever  their  nature:   electronic,  ionic  or  polyatomic) 
resist  the  change  of  the  electric  field  the  electric 
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displacement  (strain)  will  lag  the  electric  field  (stress) . 
A  lag  in  the  displacement  relative  to  the  voltage  change 
results  in  a  hysteresis  effect  and  the  scan  direction  will 
shift  the  laser  signature  in  the  direction  of  scan,  indepen- 
dent of  the  scan  rate.   From  Figures  19  -  22  nine  laser  lines 
were  obtained  in  the  laser  signature  and  the  free  spectral 
range  of  the  laser  is  approximately  1.4  kv  which  agrees  with 
the  value  found  in  Reference  13. 

Looking  at  the  shape  of  the  laser  signature  it  is 
apparent  that  none  of  the  lines  are  allowed  to  lase  over 
their  full  width  before  the  neighboring  line  begins  to 
oscillate.   The  mode  competition  or  the  fact  that  the 
oscillations  at  various  transition  frequencies  compete  for 
the  available  gain  are  quite  prevalent  in  a  CO   laser  and 
help  to  insure  single  wavelength  operation.   Throughout 
the  experiment  only  one  vibrational-rotational  line  was 
observed  to  lase  at  a  time.   The  laser  oscillated  in  the 
lowest  order  transverse  mode  TEM   throughout  the  full  range 
of  the  PS  voltage. 

The  method  of  scanning  through  the  laser  signature  can 
be  used  in  the  stabilization  system  shown  in  Figure  14  to 
find  the  required  P(20)  transition  prior  to  initiating  the 
feedback  control. 
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C.   STARK  CELL  OPERATION  AND  TEST 
1.   Stark  Cell  Design 

The  Stark  cell  described  in  Reference  13  was 
expected  to  operate  satisfactorily  in  the  staDiiization 
system-   The  cell  was  designed  to  eliminate  a  persistent 
problem  of  breakdown  in  previous  cells.   All  metal  surfaces 
were  isolated  from  the  low  pressure  gas  filled  channel 
through  which  the  laser  beam  would  pass.   The  cell  was  made 
to  accomodate  one  inch  NaCl  windows  at  each  end  and  had  a 
single  port  used  to  fill  the  cell  with  the  NH  D  gas. 
Fabrication  of  the  cell  resulted  in  an  uneven  plate  separa- 
tion of  approximately  7mm  with  1mm  of  plexiglass  between 
each  plate  and  the  channel.   The  cell  was  tested  to  3000 
volts  with  no  breakdown  and  held  a  good  vacuum.   Figure  2  3 
shows  the  cell  design  with  pertinent  specifications. 

As  discussed  in  Section  II.  D.  3  the  rest  frequency 

of  the  (0,  4    IM  =  41)  ►  (1,  5    IM  =  51)  transition  of 

NH  D  is  shifted  approximately  1700  MHz  into  coincidence 
with  the  10.5 /^m  P(20)  line  by  applying  a  Stark  field  of 
3570  ^/cm.   The  voltage  required  to  produce  the  Stark  field 
is  given  by 


V  =  -  J  £  •  dt  (49) 


where  V  is  the  difference   in  potential  between  the  plates, 
E  is  the  electric  field  strength  between  the  plates  and  d 
represents  the  path  of  integration  taken  between  the  plates 
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Although  the  simplest  path  between  the  plates  is  a 
perpendicular  straight  line  the  electrostatic  force  is 
path  independent  since  each  plate  is  an  equipotential 
surface.   Applying  (49)  to  Figure  24  for  a  straight  perpen- 
dicular path  from  the  lower  plate  (L)  to  the  upper  plate  (u) 
yields. 


V 


^U  ^U  Q  y-U  -. 

=  -)^E    •    d      =-l^E   cos    180      d      =    (_E   dje   =   E    [d- (a+b)J+E[a+bl 

L  '  L  ^  Li  O  ^ 


\fhere  E   is  the  electric  field  with  no  dielectric,  d  is  the 
o 

separation  between  plates,  a  and  b  are  the  dielectric 

thickness  at  the  plates  and  E  is  the  electric  field  inside 

the  dielectric.   Since  E  can  be  related  to  E  by  the  relation 

o 

E^ 
E  =   °/K  (51) 


where  K  is  called  the  dielectric  constant,  (50)  can  be 
written  as 

V  =  E^  [d  -  (a+b)]  +  E^  (  ^  )  (52) 

Substituting  in  the  values  for  the  Stark  cell  shown  in 
Figure  21  and  letting  E  =  3570   /cm.   The  approximate  Stark 
cell  voltage  to  bring  one  component  of  the  center  absorption 
frequency  of  NH-D  into  resonance  with  the  P(20)  line  is 
V  =  2030  volts.   Here  d  =  7mm,  a  =  b  =  1mm  and  K  =    2.7  -  3.2 
for  plexiglass  (methyl  methacrylate) .   Although  this  is  a 
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relatively  high  applied  D.C.  voltage,  it  is  within  the 
breakdown  limits  specified  for  the  cell. 


PLATE  CONFIGURATION  IN  STARK  CELL 
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2 .   Stark  Cell  Absorption  Lineshape 

Experimentation  to  determine  the  absorption  lineshape 
and  modulator  performance  as  a  function  of  pressure  began  by 
measuring  the  relative  intensity  out  of  the  Stark  cell  as  a 
function  of  pressure  with  no  electric  field  across  the  plates, 
The  attenuation  of  monochromatic  radiation  follows  Beer's  law 


1^  I(V)  =  -X(V)  I(V)  (53) 


from  which 


-  MX 

I   =  I   e  54) 

X     o 


where  x<.  is  a  linear  extinction  coefficient  which  is  the  sum 
of  the  coefficients  of  total  absorption  and  total  nonforward 
scattering.   The  relative  value  of  /^   depends  strongly  on  the 
density  and  molecular  composition  of  the  gas.   Absorption  is 
of  greatest  importance  in  the  infrared  (particularly  in  the 
8  -  14 //m  wavelength)  .   Rayleigh  scattering  because  of  its 
strong  dependence  on  wavelength  should  be  unimportant  for 
wavelengths  longer  than  the  visible  range  and  no  Mie  scattering 
should  be  present.   Since  the  density  is  proportional  to  the 
pressure  of  the  gas  a^   is  dependent  on  the  pressure. 

The  experimental  arrangement  is  schematically 
represented  in  Figure  25.   The  laser  was  tuned  and  maintained 
at  the  center  of  the  P(20)  transition  by  keeping  the  power 
meter  at  its  maximum  value  and  monitoring  the  spectrum 
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analyzer  for  the  P(20)  line.   The  cell  was  evacuated  to 

2  X  10    torr  and  used  as  a  zero  pressure  reference.   The 

Stark  cell  output  intensity  was  measured  by  the  pyroelectric 

detector.   The  y  axis  of  the  x  -  y  plotter  sampled  the 

relative  intensity  of  the  chopped  Stark  cell  output  via  the 

detector,  amplifier  and  filtered  rectifier.   The  x-axis  of 

the  plotter  was  manually  set  to  a  voltage  proportional  to 

the  pressure  of  the  NH  D  gas  in  the  cell  with  the  switch 

2 

closed  to  position  (a) .   NH^D  gas  was  obtained  by  mixing 
a  1  to  1  ratio  of  NH  and  a  commercially  prepared  sample  of 
ND^  in  the  cell.   The  gas  was  allowed  to  reach  an  equilibrium 
condition  before  the  relative  intensity  measurements  were 
taken.   Figure  26  shows  the  experimental  relative  intensity 
of  the  P(20)  transition  vs  NH^D  pressure. 

An  X  -  y  plot  of  the  relative  modulated  intensity 
vs  the  applied  Stark  cell  bias  voltage  should  provide  the 
Stark  shifted  absorption  lineshape.   The  experimental 
arrangement  is  the  same  as  Figure  2  3  except  that  the  D.C. 
Stark  cell  voltage  is  voltage-divided  and  fed  to  the  x  axis 
of  the  plotter  with  the  switch  in  position  (b) .   The  laser 
is  set  to  oscillate  on  the  P(20)  line  and  monitored.   The 
NH2D  in  the  Stark  cell  is  mixed  in  a  1:1  ratio  as  described 
above.   The  applied  voltage  to  the  cell  was  manually  scanned 
from  0  -  3000  volts.   Figures  27  -  34  show  the  experimental 
relative  intensity  vs  applied  Stark  cell  voltage  for  several 
NH^D  pressures  and  scan  directions.   An  unexpected  discontin- 
uity evident  on  the  graphs  and  while  monitoring  the  pyroelectric 
detector  output  directly  on  an  oscilloscope  showed  up  at 
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RELATIVE  INTENSITY  vs  NH2D  PRESSURE  W   STARK  CELL 
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pressures  above  5  torr.   This  abrupt  change  in  the  modulation 
of  the  laser  made  frequency  stabilization  impossible. 
Repetition  of  the  experiment  resulted  in  similar  curves 
with  no  significant  differences.   It  is  evident  from  the 
curves  in  Figures  29  -  34  that  the  scan  direction  signifi- 
cantly effected  the  shape  of  the  curve  and  location  of  the 
discontinuity,  that  the  discontinuity  occurred  at  different 
applied  voltages  to  the  Stark  cell  for  different  NH-D 
pressures  in  the  cell  and  that  a  substantial  amount  of 
modulation  of  the  beam  (as  much  as  50%)  could  be  achieved 

using  NH  D  and  the  Stark  effect. 
2 

Closer  inspection  of  Figures  29  -  34  indicated  that 
the  curve  began  to  repeat  itself  after  a  discontinuity 
occurred  independent  of  the  scan  direction.   Figure  35  is 
an  experimental  plot  of  the  relative  intensity  of  the  P(20) 
laser  transition  vs  voltage  applied  to  the  Stark  cell  with 
the  curves  for  both  scan  directions  (a  -*•  0  -  3000v  and 
b  — *-  3000  -  Ov)  on  the  same  graph.   It  shows  clearly  the 
effect  of  scan  direction  on  the  plot  and  repetition  of  the 
curve  shape  after  the  discontinuity. 

This  repetition  effect  was  studied  by  reversing  the 
scan  direction  of  the  voltage  applied  to  the  Stark  cell 
after  each  discontinuity  yielding  the  plots  shown  in  Figures 
36  and  37.   The  curves  did  repeat  themselves  with  basically 
the  same  shape  and  with  the  discontinuities  occurring  at 
different  Stark  cell  voltages  and  relative  intensities.   The 
curves  of  Figures  36  and  37  could  be  retraced  for  both 
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directions  of  scan  provided  the  voltages  corresponding  to 
the  discontinuity  were  not  reached.   Thus  the  shape  of  the 
irradiance  is  independent  of  scan  direction. 

In  order  to  observe  the  operation  of  the  dither 
mechanism  an  A.C.  voltage  was  applied  to  one  plate  while  the 
D.C.  voltage  on  the  cell  was  maintained  at  a  voltage  where 
sufficient  modulation  was  observed  but  below  the  point 
where  the  discontinuity  occurred.   A  200Hz  A.C.  signal  was 
supplied  to  one  of  the  Stark  plates  from  an  audio  oscillator 
amplified  to  60  volts  vms .   The  laser  was  set  to  oscillate 
on  the  P(20)  transition  and  the  chopper  was  removed  from 
the  experimental  arrangement  in  Figure  25.   The  signal  from 
the  pyroelectric  detector  was  fed  directly  into  the 
oscilloscope.   A  200Hz  sinusoidal  signal  was  observed.   This 
200Hz  sinusoidal  signal  was  lost  when  the  laser  was  tuned 
off  the  P(20)  line  by  changing  the  applied  voltage  to  the 
piezoelectric  driver  in  the  laser.   A  scan  through  the 
laser's  spectral  signature  revealed  the  expected  results 
of  modulation  of  the  P(20)  and  P  (14)  laser  lines  by  NH-D. 
A  200Hz  sinusoidal  signal  was  evident  on  the  oscilloscope 
for  both  these  transitions  with  the  amplitude  of  the  signal 
for  the  P(14)  line  being  about  50%  of  the  amplitude  of  the 
signal  for  the  P(20)  line  for  the  preset  Stark  cell  voltage. 
Thus  the  Stark  cell  was  actually  modulating  the  laser  beam 
and  could  produce  a  suitable  signal  for  phase  detection  if 
the  discontinuity  could  be  corrected. 

Attempts  to  determine  the  cause  of  the  discontinuity 
were  unsuccessful.   The  voltage  across  the  cell  was  equal 
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to  the  supply  voltage  (0  -  3000VDC)  with  no  detectable 
voltage  drop  existing  that  might  have  caused  a  discontinuity. 
No  current  could  be  detected  through  the  cell  which  would 
have  indicated  ionization  of  the  gas  or  cell  breakdown.   In 
order  to  determine  if  perhaps  the  beam  was  being  steered 
off  the  detector  by  a  refractive  index  gradient  in  the  cell, 
rather  than  the  gas  modulating  the  beam,  a  detector  was 
scanned  across  the  output  of  the  cell  when  the  applied 
voltage  on  the  cell  was  at  a  point  where  it  was  able  to 
modulate  the  P(20)  line.   The  pyroelectric  detector  was 
mounted  on  a  moveable  table  and  scanned  across  the  Stark 
cell  output  signal  giving  a  relative  plot  of  intensity  vs 
position.   The  signal  from  the  detector  was  amplified, 
filtered  and  rectified  and  then  fed  into  the  y  channel  of 
a  time  base  plotter.   The  scan  rate  of  the  detector  was 
fixed  so  that  the  time  base  could  be  translated  into  a 
position  measurement.   Figure  38(a)  shows  the  beam  shape 
out  of  the  cell  as  a  function  of  position  for  zero  voltage 
applied  to  the  cell  and  Figure  38 (b)  the  beam  shape  out  of 
the  cell  when  the  beam  was  being  modulated  with  an  applied 
Stark  cell  voltage.   The  intensity  of  the  beam  is  reduced 
with  no  noticable  shift  in  position  of  the  beam.   The 
detector  was  positioned  at  the  point  of  maximum  intensity 
and  a  D.C.  bias  voltage  was  applied  to  modulate  the  base. 
The  detector  then  moved  to  determine  if  the  beam  was  being 
steered  to  the  left  or  right,  up  or  down.   Figure  39  is  a 
graph  of  intensity  vs  time  and  shows  that  the  beam  was  not 
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EXPERIMENTAL  PLOT  TO  DETERxMINE  IF  A  REERilGTIVE 
INDEX  GRADIENT  EXISTED  IN  THE  STARK  CELL, 
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being  steered.   The  intensity  as  the  detector  was  moved  did 
not  exceed  the  modulated  intensity.   Again  the  discontinuity 
is  evident  as  the  applied  voltage  was  increased. 

The  movable  detector  was  placed  at  various  positions 
in  Figure  2  2  to  determine  if  the  beam  shape  was  due  to 
reflections  in  the  cell  or  prior  to  the  cell  itself.   Figure  40 
shows  the  relative  intensity  beam  profile  of  the  laser  at 
three  positions  prior  to  the  cell:   40(a)  directly  out  of 
the  laser  prior  to  beam  splitter  No.  1,  40(b)  between  beam 
splitters  No.  1  and  No.  2  and  40(c)  between  beam  splitter 
No.  2  and  the  cell.   It  appears  that  the  germanium  beam 
splitters  cause  some  distortion  of  the  beam  and  contribute 
to  the  distorted  shape  of  the  beam  out  of  the  Stark  cell. 
The  profile  to  the  left  of  the  main  profile  in  Figures  38(a) 
and  38(b)  was  probably  caused  by  reflection  from  back  surface 
of  the  parallel  sided  Ge  beam  splitter  since  it  appears  to 
be  a  mirror  image  of  the  main  profile. 
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3 .   New  Stark  Cell  Design 

In  order  to  determine  if  the  discontinuity  was  due  to 
the  nonuniformity  of  Stark  cell  plate  spacing  in  the  cell  or 
the  cell  design  itself,  another  Stark  cell  patterned  after 
one  described  in  Reference  31  was  built.   The  new  Stark  cell 
design  with  its  specifications  is  shown  in  Figure  41.   Its 
design  was  an  attempt  to  reduce  the  applied  voltage  require- 
ments by  reducing  the  plate  spacing  to  3mm  and  eliminating 
the  need  for  spacers  to  keep  the  plates  parallel  by  attaching 
the  plates  to  plexiglas  holders  with  screws .   This  design 
allows  for  removal  of  the  plates  for  modification,  cleaning  or 
substitution  of  a  new  design  for  the  plates.   The  plate  thick- 
ness can  be  used  to  vary  their  separation  in  the  cell.   The 
edges  and  ends  of  the  plates  were  rounded  to  help  prevent 
breakdown.   The  end  caps  hold  the  NaCl  windows  and  simplify 
holding  of  the  windows  when  modifying  the  cell.   The  cell  held 
a  good  vacuum.   The  cell  was  filled  with  approximately  2  torr 
of  NH  D  with  breakdown  and  ionization  of  the  gas  occurring  at 
1000  volts.   The  breakdown  was  observed  to  occur  along  the 
edges  and  especially  at  the  ends  of  the  plates.   It  is  believed 
that  the  inability  to  machine  round  edges  is  the  primary  cause 
for  breakdown.   Investigations  into  fabrication  of  plates  with 
round  edges,  insulating  the  edges,  use  of  geometries  other  than 
strictly  flat  plates  and  electroplating  of  the  plates  on  an 
insulator  need  to  be  conducted  to  produce  a  viable  cell  for 
frequency  stabilization  using  Stark  cell  techniques.   The  pre- 
sent cell  can  be  used  to  test  plate  designs  and  fabrication 
techniques  without  having  to  build  a  new  cell  each  time. 
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NEW  STAM  CELL  DESIGN  WITH  INTERCHAI^GEABLE  PLATES 


NaCl 
WINDOW 


VOLTAGE 

GOMSCTOR 


PARALLEL 
PLATES 


GAS 

INLET 


T 

5  cm 

1 


o        o 

^  ^•'  A 

O     ;  ;     O 

- 

J 

1 

T 


3  mm 


SPECIFICATIO 

N3 

PLATE: 

MATERIAL 

ALUI'IINUM 

LENGTH 

15.0  cm 

WIDTH 

2,5  cm 

SEPARATION 

3.0  mm 

CELL: 

MATERIAL 

PLEXIGLAS 

FIGURE  41 


101 


IV.   CONCLUSIONS 

A  Stable  high  voltage  power  supply  (Sylvania  941S)  run 
from  a  stable  line  regulator  was  found  suitable  as  a  power 
source  for  the  Sylvania  CO  Model  941  Laser  in  a  Stark  cell 
stabilization  system. 

A  servo  loop  integrator  was  designed  and  built  allowing 
the  correction  signal  to  the  piezoelectric  driver  to  be 
proportional  to  the  time  integral  of  the  error  signal  pro- 
duced by  phase  sensitive  detection.   The  integrator  used  in 
the  feedback  system  should  provide  increased  stability. 

A  method  for  scanning  through  the  laser's  spectral  signa- 
ture to  find  the  required  P(20)  transition  incorporating  the 
same  integrator  used  in  the  feedback  stabilization  system 
was  demonstrated  and  discussed. 

A  Stark  cell  was  tested  to  determine  its  characteristics 
for  use  in  the  stabilization  system.   It  was  found  that  a 
discontinuity  existed  in  the  cell's  lineshape  function  which 
made  it  unsuitable  for  use  in  a  phase  detection  stabilization 
scheme.   The  cause  of  the  discontinuity  could  not  be 
determined.   A  voltage  drop  across  the  cell,  ionization  of 
the  NH2D  gas,  and  beam  steering  caused  by  an  index  gradient 
being  created  in  the  cell  were  eliminated  as  possible  causes. 

A  second  new  Stark  cell  with  the  capability  of  modifying 
or  changing  the  cell  plates  was  designed  and  built  to  deter- 
mine if  the  discontinuity  found  in  the  first  cell  was  inherent 
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in  the  cell  design  or  caused  by  a  factor  outside  the  cell. 
Breakdown  occurred  in  the  second  cell  prior  to  the  field 
strength  required  for  modulation  making  the  second  cell 
unsuitable  for  use  in  the  system.   The  fabrication  of  the 
plates  in  the  new  cell  was  the  major  cause  for  breakdown. 
Further  investigation  needs  to  be  conducted  in  the  Stark 
cell  plate  design  and  fabrication  to  prevent  breakdown 
inside  the  cell. 

Production  of  an  error  signal  using  a  pyroelectric 
detector  was  demonstrated  using  the  first  Stark  cell  but 
was  unsuitable  for  phase  sensitive  detection  due  to  the 
discontinuity  in  the  lineshape.   It  is  expected  that  the 
successful  fabrication  and  test  of  a  Stark  cell  will  produce 
a  Stark  stabilized  CO   laser  in  the  assembled  system.   This 
system  frequency  stabilized  on  the  P(20)  transition  is 
expected  to  be  operational  in  both  the  laboratory  and  the 
field. 
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APPENDIX  A 
MODEL  941  SYLVANIA  L?\SER  SPECIFICATIONS 

Output  v/avelength:   10.6 /cm 

Power  output:   3  watts 

Beam  divergence:   4  mrad 

Beam  diameter  to  l/e  point:   4  mm 

Beam  polarization:   vertical 

Cavity  length:   45  cm 

Bore  diameter:     6  mm 

Cooling  requirements:   tap  water  @  15  to  30  ml/s 

Output  mirror:   70%  reflecting,  dielectric  flat 

Piezoelectric  mirror:   100%  reflecting,  300  cm  radius  of 

curvature 

Piezoelectric  ^   parameter  as  given  by 

_3 
Sylvania:   -4  x  10     m/volt 

Laser  voltage:   14  kv  (3  7.6  ma 

Tube  voltage:   9.5  kV  @  7.6  ma 

Series  ballast  resistance:   600  kO. 

Equivalent  tube  resistance  @  7.6  ma:   1.27  M 

Pressure  broadened  linewidth:   100  MHz 

Doppler  broadened  linewidth:   50  MHz 

^   mode  separation:   333  MHz 
2L 

Cavity  linewidth:   19  MHz 
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APPENDIX  B 


OPTICAL  ENGINEERING  CO   LASER  SPECTRUM  ANALYZER 

2 


The  CO   Laser  Spectrum  Analyzer  is  an  infrared  grating 
spectroscope  designed  simultaneously  to  monitor  and  visually 
display  through  the  use  of  a  calibrated  thermal  sensitive 
screen  which  darkens  at  points  corresponding  to  the  wave- 
length at  which  the  laser  is  operating.   The  optical  design 
of  the  instrument  is  basically  an  Ebert  Mount  using  a  plane 
grating  and  a  spherical  mirror.   Folding  mirrors  are  used  to 
reduce  the  length  and  the  design  modified  so  that  the  display 
screen  can  be  located  near  the  top  of  the  instrument. 
Incorrect  positioning  of  the  instrument  so  that  the  laser 
light  does  not  enter  within  the  acceptance  cone,  will  not 
create  incorrect  readings  but  simply  cause  the  light  to 
miss  some  element  in  the  optical  path  and  not  reach  the 
display  screen. 

Figure  IB  is  a  side-view  of  the  spectrum  analyzer 
showing  its  essential  components.   In  normal  operation 
the  portion  of  the  laser  beam  which  passes  through  the 
slit  will  strike  a  flat  mirror  which  folds  the  optical 
beam  and  introduces  a  slight  vertical  angular  component  to 
the  direction  of  the  beam.   The  light  then  strikes  a  concave 
mirror  which  collimates  the  divergence  introduced  by  the 
diffraction  of  the  slit.   The  light  then  is  reflected  to 
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the  diffraction  grating.   The  grating  is  rigidly  mounted 
on  the  flat  mirror  and  held  at  a  fixed  angle.   Only  the 
light  which  is  diffracted  into  the  first  order  of  the 
diffraction  grating  continues  on  to  strike  the  curved 
mirror  for  a  second  time  which  focuses  the  image  on  the 
display  screen.   Each  laser  wavelength  leaves  the  grating 
at  a  unique  angle  and  is  imaged  at  a  different  point  on  the 
calibrated  scale. 

The  specifications  for  the  spectrum  analyzer  are  given 
in  Table  I-B. 
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TABLE  I-B 
SPECIFICATIONS 
The  following  is  a  list  of  the  specifications  for  the 
Model  16A  CO  Laser  Spectrum  Analyzer. 

Wavelength  range:   9.1  to  11.3 

Resolution:   .003  (0.3  cm  ) 

Maximum  power:   75  watts 

Minimum  power  (focused):   0.01  watts 

Screen  responce  time:  h   second 

Grating  efficiency:   90%  horizontal  polarization 

30%  vertical  polarization 

Optical  design:   3/4  meter  modified  Ebert 

Total  optical  path  length:   2.3  meters 

o  o 

Acceptance  angle:   2.5   horizontal,  1.5   vertical 

Slit  width:   adjustable  -  closed  to  0.6  mm 

Weight:   12  lbs. 

Length:   18" 

Width:   9.9" 

Minimum  height:   6.5" 

Elevation  of  slit:   2.5"  to  3.8"  (adjustable) 

Electrical  power:   115  V,  60  Hz,  10  watts 

12  16 

C  O^  vibrational-rotational  lines    ^      q 

identified:   00  1  -  02  0 

P2  to  P62  and  R2  to  R60 

00°1  -  10°0 

P2  to  P56  and  R2  to  R60 

Ol-'-l  -  11^0 
P19  to  P45 
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APPENDIX  C 
SERVO  LOOP  INTEGRATOR 
The  output  of  an  ideal  integrator  is  a  voltage  pro- 
portional to  the  integral  of  the  input  voltage;  or  in  other 
words,  the  output  is  proportional  to  the  product  of  the 
amplitude  and  duration  of  the  input.   The  integrator  performs 
this  mathematical  operation  on  an  instantaneous  basis.   The 
result  is  an  output  exactly  proportional  to  the  area  under 
a  wave  form. 

The  circuit  in  Figure  B-1  performs  integration  using  a 
yuA  747  dual  frequency  compensated  operational  amplifier  to 
force  the  same  current  through  both  R  and  C  .   The  voltage 
across  the  feedback  capacitor  is  related  to  capacitor  current  by 


v^ 37|if'3t         (ic) 


Since  the  circuit  causes  i  to  equal  the  input  current 

V. 


(i,  =    Vr^) 


^1  Cf  J   ""i 


V  =  V  =  -  -r — :: —  \  V.  dt  (2c) 

c    o      R- 


The  gain  of  the  circuit  is  given  by  /R  C^.   Thus  the  output 
voltage  will  change  by  -  /C  R.       /sec  for  each  volt  of  output 

The  bias  voltage,  V,  ,  and  variable  voltage  offset  null 
resistor,  R  ,  the  drift  inherent  in  the  integrator.   The 
switch,  S  ,  will  reset  the  integrator  to  zero. 
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SERVO  LOOP  EITEGRATOR  CIRCUIT 


SPECIFICATIONS: 


OPERATIOML  AliPLIFIER:     >f  A  747  DUAL  FREQUENCY  COMPENSATED 


R^   .   Rp  =  R^  =   1  MA 

R^  =   10  kil 

R^  =   1   k/l 

C^  =    1  >u? 

R     s   10  kA  (variable) 

Gain  =  "  /R^C^  =  -1 

V.  =  input  voltage 

V  =  output  voltage 

V,  =  bias  voltage 

S  =  reset  switch 
r 


FIGURE  C-1 
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